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Abstract 

Recently there has been a growing interest in Glycol fluid solubility 
in water. The applications of Glycol- Water fluids is in the biomolecu- 
lar field, physical chemistry, and engineering, from protein dynamics 
to power production and storage devices technology. In this paper a 
theoretical model for nanometer scale or mesoscale interacting fluids 
is derived for the interacting Glycol and water molecules. The theory 
is a semi-classical interacting potential theory which utilizes effective 
potentials such as of the Lennard-Jones type. The statistical distri- 
butions of the fluids are derived from information theory and tem- 
perature dependent distributions are obtained for the Lennard-Jones 
potential interactions of the fluids. 

1 Introduction 

The applications of Glycol- Water fluids is in the biomolecular field, physical 
chemistry, and engineering, from protein dynamics to power production and 
storage devices technology. In this paper a theoretical model for nanometer 
scale or mesoscale interacting fluids is derived for the interacting Glycol and 
water molecules. The theory is a semi-classical interacting potential theory 
which utilizes effective potentials such as of the Lennard-Jones type. The 
statistical distributions of the fluids are derived from information theory 
and temperature dependent distributions are obtained for the Lennard-Jones 
potential interactions of the fluids. 

The full description of the fluids of water and Glycol would require a fully 
quantum mechanical theory. Such a theory exists and will be briefly discussed 
in the derivation. However for practical applications, simplifications of the 
computation required to model the biomolecular dynamics of Glycol- Water 



are sought. Simplifications are possible and are motivated by the success of 
Lennard- Jones molecular dynamics of fluids and biomolecular interactions, 
at nanometer scales. Such models the utilize Lennard- Jones type of poten- 
tial interactions for example are incorporated in computational molecular 
dynamics such as the NAmd Nanomolecular dynamics source code free dis- 
tributions of the University of Illinois molecular dynamics of biomolecules 
project [2]. 

The theoretical model will be derived for the two fluids as independent 
fluids. Then the interacting binary fluids model will be developed and dis- 
cussed. The resulting statistics will be numerically simulated and compared 
with experimental results published by researchers that have conducted de- 
tailed experimental studies of Glycol- Water interactions. 

2 Interaction between differing molecular species 
and phases 

The Glycol interacts with water by Hydrogen- Oxygen bonding. The Water 
is attracted to the Oxygen atoms and Hydrogen bonds are formed in the case 
of ethylene oxide adducts dissolved in water. 

-CH 2 -O- CH 2 - 
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A 1980 study by steuter, Mozafar and Goodin [13] also has shown that some 
of the ethylene adducts in aqueous solutions 'may be present as cation-active 
polyoxonium compounds' 

* (2) 
[-CH 2 -0-CH 2 -}+(OH)- 

With these two types of reactions (1) and (2) co-existing in an equilibrium 
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3 Water model 



Water is a fluid that is the universal solvent of nature . It is a fluid at 
standard temperature and pressure. The fluidity of water is of a complex 
arrangement of dipole-dipole interactions of the polar molecule. Water has 
a triangular shape with a bend angle of a — 104.5°. Water is comprised of 
H — O — H connected by covalent bonds and forms Hydrogen bonds readily 
with most molecular material. 

Water can be described by a Hamiltonian for fluids. The Hamiltonian is 
composed of the Hydrogen Hamiltonian and the Oxygen Hamiltonian and the 
interacting Hydrogen-Oxygen potential . The Hydrogen potential is Hu = 
h e - + h p + + h e - e - + h p +- p + — h e - p + 

Ho = J2 e k ct k c k + Y l k'° ] k' c k' + Y t m, n)c\c^rhc n (3) 

k k' k,l,m,n 

The gaseous atoms due to attraction form molecular bonds and become H 2 . 
To illustrate the more complex interactions and the method of description 
we discuss the interaction formalism and the phase transition theory, which 
we will apply to the more complex H 2 in a similar approach, and to the 
Glycol molecular fluid and the Glycol- Water binary fluid. ' 
The Hamiltonian for Hydrogen can be written as a matrix 

jj | h e - e - h e -p+ 

\hK- P + h p +_ p + 

where we have partitioned the Hamiltonian, a freedom we have in this and in 
all subsequent work, to include the free electron and electron-electron inter- 
action as /i e -_ e -, similarly for the free proton and proton-proton interaction, 
and the interactions between the electron-proton as the off diagonal inter- 
action terms. Note that the electron-proton interaction is of negative sign, 
a potential that lowers the energy to a more stable configuration, and also 
that the off diagonal terms are hermitian conjugates and transposes of each 
other. 

The Hydrogen molecule is then described by the long range order param- 
eters, the diagonal DLRO parameter useful for observables such as the num- 
ber of electrons and protons, as 7V( e -w p +) = iV e - + N p +, and the off-diagonal 
ODLRO useful for interacting formation of Hydrogen and Hydrogen Molec- 
ular gas as N h +h 2 — N e -_ p +. These are obtained from the retarded and 
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advanced r, a Green's functions of the Hamiltonians, which are for example 
for the electrons q r t _ = =p- ; , — and n is an infinitesimal that 

shifts the singular denominator such that outgoing (retarded) wave solutions 
for example are obtained. Also, the Fermi-Dirac distribution f(oo) is a tem- 
perature dependent distribution that is a step function at zero temperature 
and is a limit on the energy integration N = i J g r (k,u)dkdu by the disper- 
sion relation. Other observables, such as energy < E > and rates of reaction 
can also be obtained from this matrix formulation and the ODLRO as the 
derivative of the time dependent observables as djV ^ g2 = J^i J g r (k,t — t')dk. 

The phase transition of the molecular Hydrogen is described by another 
matrix interaction, the interaction being of the fluid and the lattice phases 
for example. Writing the Green's functions matrix as 

hi- J 9f-S 9f-g\ 
^lattice = 9s- f 9S-S 9S-g (5) 
\9g-f 9g-S 9g-gJ 

The Hamiltonians are simplified Hydrogenic Hamiltonians h g ^ g = E e (^) c ^ c fc+ 

k 

E <f>(k,l,m,n)c^cic'rhCfi, and the lattice Hamiltonian is shifted by con- 

k,l,rn,n 

stant wave vectors corresponding to the lattice periodicity, geometry and 
symmetry as h S -s = E e (^) ct fc+Q c fc+^ _ E _ 0(£ h™., ^ k+Qo !^ m-Q' c n- 

k+Qo k+Q ,l,m—Q' ,n 

The fluid phase Hamiltonian is similar as the lattice Hamiltonian however 

— # — # 

the constancy of the lattice vector is relaxed to Q — > Q a variable. Other 
simplifications are possible, such as solid phase tight binding Hamiltonians, 
metallic effective mass Hamiltonians. 

The formation of water is a similar interaction description from this for- 
mulation of many-particle physics of interacting atoms and molecules and 
phase transitions . The formation of the water molecule in the fluid phase 
is approached from the gaseous or vapor phase for definiteness. The Hydro- 
gen gas of Hydrogen atoms and H 2 molecules is for example allowed into an 
Oxygen environment and the two gases interact as due to Oxygen's effective 
electron affinity. The matrix Green's function for the formation of Water 
H 2 vapor then is written as 

f 9h-h gH-o\ 
9(H 2 o) vapor - [ t „ J 
\9o-h 9o-o) 

and the Hamiltonians are written as discussed with the appropriate symme- 
try, approximation as needed to the required accuracy desired. 



(6) 
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The fluid phase is obtained as before , with the Hamiltonians including 
interaction terms of the phase symmetry, 



_ / 9v—V 9v — f \ frj\ 

9{H 2 0) vapor - fi uid — I f I \l) 

\9 f-v 9f~fj 

The Hamiltonians are of interactions between the vapor-fluid phases, 
and the interactions are of the gas potential with the fluid potential vec- 
tors with the 'quasi-periodic' lattice vectors shifted by the variable vector 
< $ >= (f>(k,Q,l,k',l')c^ +( ^cpaiCp and here this denotes that a 

k+Q,l,k'X' 

vapor 'particle' is created at I and propagated through the system , interact- 
ing with the fluid 'particle' of k + Q the two particles then being removed 
at their final momentums of k',l' respectively and here we have written the 
fluid symmetry wave vector explicitly. 



4 Lennard-Jones theory of multi-phase flu- 
ids. 

The description of water from the fully quantum mechanical can be simpli- 
fied. An indication of this was already given when we discussed the level 
of accuracy, approximations utilized to achieve this accuracy , such as the 
use of Hydrogenic wave functions, tight binding Hamiltonians, effective mass 
Hamiltonians... An approximation that can be made is an approximation of 
the dynamics as quantum effects become screened. As the number of atoms 
and molecules increases, the screening that occurs makes it less necessary to 
describe the interactions at the fully quantum mechanical description. The 
description of the interactions can be made then as semi-classical interac- 
tions. The most well-known of theses descriptions being the Lennard-Jones 
potential and its generalization 0(r) = fji — ^ description of interactions, 
which is a screened Coulomb potential with the Pauli exclusion principle 
accounted for by the infinite amount of energy needed for two particles to 
occupy the same space, this being a mathematical description of repulsion 
of a sort. The Lennard-Jones potential is further simplified by the approx- 
imation by a harmonic oscillator potential in the vicinity of the minima of 
the potential well. However this approximation we consider as not being of 
sufficient accuracy to describe the dynamics of the fluid phase of water, and 
the subsequent description of the dynamics of Glycol. 

Nanofluid droplets and the physics of fluids at the nanometer scale and 
at the meso scale between the fully quantum mechanical and the classical 
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continuum of hydro-dynamics have been described reasonably successfully in 
the literature by semi classical methods [6], [8] . The approaches have been to 
classically treat the nanofluid droplet by approximate hydrodynamic means, 
the vapor by Boltzmann equation, and to evolve the number of particles and 
obtain rates of change in number of particles. Approximations are made such 
as Boltzmann equation [8] collision terms of a certain order, an example is 
the force collision terms up to a third order expansion by [7]. 

In hydrodynamics approaches, full analytic solutions are not obtained to 
the hydrodynamic equations... approximations of linearization and similar 
assumptions for solution are made such as in [6] , examples of approximations 
are of viscosity and compressibility. 

Another approach between the hydrodynamic approximation and the 
quantum mechanical dynamics of the mesoscale is possible. Given the success 
of recent molecular dynamics simulations in describing experimental results, 
approximating the nanofluid as a Lennard- Jones fluid, this approach seems 
reasonable. Then a similar theory of multi-phase fluid-vapor can be made 
at the Lennard- Jones approximation level that is applicable in between the 
range of applicability of the Boltzmann and hydrodynamics equations' scale 
of applicability and the fully quantum mechanical dynamics. That is, from 
the hundreds of atoms to the micrometer scale. From an analogy to the 
fully quantum mechanical phase transition theory discussed in the previous 
section, a theory of coexisting phases of nanofluid droplets can be obtained. 
The advantages would be the simultaneous description of co-existing phases 
by analogy to the matrix formulation discussed, long range order parame- 
ter approach as in the diagonal trace and summation to obtain particle or 
number of atoms in the individual phases. 

A theoretical advantage of the full quantum theory aside from increased 
accuracy at the 'cost' of increased computation time is the generality of the 
theory and the well known Hamiltonian superposition principle that allows 
us to include external perturbations and interaction terms and which allows 
us to write a theory of phase transitions as a superposition of interactions 
between the co-existing phases. From the point of view of the approach of 
Lennard- Jones fluids coexisting phases theory, the superposition principle 
is replaced by the Legendre transform of the maximum entropy method of 
theory derivation, and equivalently the information theoretic method. 

Alternatively a Lennard- Jones semi classical Lagrangian or Hamiltonian 
may be used, and dynamics with or without a heat bath computed. The 
heat bath dynamics implying fluctuations and therefore stochastic dynam- 
ics, equivalent descriptions made by Fokker-Planck equations can be ob- 
tained. These are analogous to the quantum mechanical Green's functions 
Schroedinger-like equations and can therefore be written in matrix form for 
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the co-existing phases. We derive our theory utilizing the entropy maximiza- 
tion method. 

4.1 Derivation 

The entropy state function or information measure to be used we choose to be 
the Gibbs-Boltzmann entropy, < S >= —c J P(f, t)lnP(f, t)dr which can be 
generalized to the nonextensive entropy of C. Tsallis [TTJ [TJ] for better accu- 
racy in future work.... also other novel effects such as the q-parametrization 
of the deviation of the PDF from the Gaussian to the power-law PDF can 
be advantageously made for simple interacting nonlinear dynamics. Note 
that nonlinearity due to physical interactions deviates the PDF from the 
Gaussian distribution towards the power-law distribution. The nonextensive 
statistics has been noted for theoretically producing this form of power-law 
distribution from a formally similar mathematics to the traditional thermo- 
dynamics and information theory. The connections are more profound and 
review articles and applications can be found at [11] . 

The maximum entropy approach is then of the maximization of the en- 
tropy measure constrained with observables of the system in a Legendre 
transform with Lagrange multipliers setting the weighted units to the en- 
tropic measure's. 

The choice of observables is one of choosing the minimum set of observ- 
ables that describe the physics of interest. These observables are described 
mathematically as functions of the other observables chosen as representative 
of the observable dynamics of the system. A discussion of this was made by 
the original author of maximum entropy theory [91 \TU\ E. T. Jaynes. 

We choose observables analogous to the quantum mechanical observables 
in the fully quantum mechanical theory. The momentum or kinetic energy 
observables of quantum mechanics correspond to the first and second mo- 
ments of the variables of position, the potential energy or interaction poten- 
tial terms are similarly included here though we choose the Lennard- Jones 
potential and not the Coulomb potential, The observable moments are: 

hi 

= y~] I (Axij- < Axij >) 2 P(Ax7j, Axfji, ,At)dAx~lj (8) 
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N 

M 2 (A%,) = E < (Af iy - < A%, >) 2 > 



/ (Afj/j/- < Afj/j/ >) 2 P(Ax ij ,Ax^ j >, ,At)dAxJj> (9) 

i'J' ^ 



AT i? 

oHM 2 (Ax k "k>) = E < (Ax k »k>- < Ax k „ k , >) 2 > 
fc",fc' 

E / (Ax k "k'- < Ax k n k , >) 2 P(Ax ij ,Ax^ j f, ,At)dAxJj> (10) 



7V H 

< Vh-h >- 2^ < (— ^-J > 



ATff , 

E / (^f - ^)^(A% Af )7) , At)dA^rfA%, (11) 



N 

< vo-o >= £ < (^r " ^) > 

No , 

E / (^T - ^)P(A% Ax iY , , A0dA£^A%, (12) 



i,3 



Nh,No 

,vh~o Vh-o 



E / " ^) P ( A ^' A %'> > AQdAxijdAxi.j, (13) 



<l/ OJI _o>= £ < («»^ - 222=2) > 
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< VH20-H20 > z 



Nh20 

E< 

k,i 



&H20-H20 VH20-H20 



)> 



kl 



< VmO-Y > z 



N H 20,N Y 

E < 

k,l> 



&H20-Y f]H20-Y 



Tkl^ 



> 



kl' 



N H' N ([Y] + (OH)-) 



< v, 



H-([Y] + (OH)~) 



> = 



E 



<( 



a H-([Y] + (OH)-) VH-([Y] + (OH)-) 



)> 



These observables of free 'particle' movement of trajectories and ther- 
mal fluctuations and electromagnetic interactions as described by the semi- 
classical potentials of a Lennard- Jones type whose semi-classical approxima- 
tion is general enough to have force constants that can be determined from 
experiment or alternative theoretical derivation. The interaction between 
particles of atoms forms water, and ions of Hydroxide and molecules of H 2 
and O2 and interacting H 2 molecules and interactions of water molecules 
with the Glycol's ether Oxygen atoms here denoted by species Y. 

The maximization of entropy is made with the Lagrange multipliers set- 
ting appropriate units for the observables as 5[< S >]+ 51(3((IhhM 2 + + 



(Ih-([y]+(oh)-) < Vh-([y]+(oh)-) >] == which derives the least biased dis- 
tribution matrix for binary interactions 



/ Ph-h 



H-([Y] + (OH)-) 



\Ph-([y]+(oh)- 



) 



(14) 



P([Y]+(OH)-)-([Y]+(OH)-)J 



We will specialize this to the simple case of fluid water solvent and dis- 
solved ether, the steps in derivation thus far being instructive as how to 
proceed. 



5 Lennard-Jones Water-Glycol interactions 

The simplified case of the water and dissolved adducts of the Glycol that is 
the ether and ethylene can be written as the previous derivation which we 
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omit. We write the least biased distributions 

Py&Xij, , At) = A (H2 o,h 2 o) e « « 



.,Af) = 



(ff,([y]+(Ofi)-)) 



-<A^. 



» 2 ]+a 



(H,[y]+(Off)-) 



'g,([y]+(OH)-) ''g.nyj+tofl 



The Lagrange multipliers are f3 = 2DAf which are the inverse of the diffusion 
coefficient and the Diffusion constant is by the Einstein relation proportional 
to the inverse temperature from thermodynamics. The normalization is the 
inverse of the partition function and can be evaluated in several methods, 
we will leave that for later as it can be normalized arbitrarily or to unity. 
The solution of this least biased distribution is dependent upon the Lagrange 
multipliers and if normalization is needed then integration over the variables 
to obtain the partition function. Some identities utilized for evaluation are 
A(t) = the normalization and partition functions , (3(t)Z(t) 2 = const. 
and -||| = — < (Axi'ji— < Arr^y >) 2 > and similarly for other Lagrange 
multipliers and observables. 
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5.1 Derivation of Solutions 

P(A£ W , At) = e - Z(M) - - . (16) 

The Lagrange multipliers are evaluated as discussed from the relation- 
ship between partition function and the observables = —M 2 (Axij) and 
similarly for the other observables. Note that the Lagrange multipliers of the 
second moments are related to the diffusion coefficients by the temperature- 
diffusion coefficient thermodynamics relation D = A, and < (xjj) 2 ) >= 2Dt. 
The analytic solution of interest here is the one of the transformation of the 
Lennard- Jones potential to a drift coefficient term. This is done by uncom- 
pleting the square of the Gaussian 

p(X) = e ( x2 ~ 2X < x >+< x > 2 )/ 2Dt (17) 

and comparing it with the dependent terms of the de-mean (or constant 
or linear and renormalized mean) Gaussian 

P(r l3 , ....) = e -r i ^ r 2n ] <n 3 >+<n ] >>)- 1 LJ(n ] )) e g ^ 

The Lennard- Jones potential terms can be seen to be terms that can be 
transformed to a drift coefficient or the first moment by un-completing the 
square. We obtain 

a(fij) =< nj >= nj + ^nf + jjLJ(rij). (19) 

5.2 Analytic distribution functions 

With this the r\j dependent terms become a simple e -/ 3 ( r »i- a ( r ij. A ')) 2 Gaus- 
sian which by the identity ^f- = -M 2 (A£^) gives the simple diffusion- 
temperature relation (3 = with the temperature from thermodynamics 
proportional to the diffusion constant. 

The resulting PDF is a Gaussian with the coordinates expressed in spher- 
ical coordinates for simplicity. The integration for the partition function or 
normalization can be accomplished, however we use the relationship (3Z(t) 2 = 
2n and the resulting distribution is 

P(r lv ....) = A r (t) 2 . (20) 

y/AirDt 

where the derivation implies 
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a) the Lagrange multiplier second moment are related as < r >= 2Dt = 

1 



b) the inverse partition function and the normalization are therefore also 
Z r (t) = yj An Df 2 

c) the partition function Z(t) = z r (t)zg(t)z^(t) can be factorized and, 

d) Lagrange multiplier and factorized partition functions are also con- 
stants as the overall identity (3Z 2 = constant. 

Note that aside from the r^ terms which have the complicated Lennard- 
Jones potential dependence and its radial coordinates dependence Eq.(19) 
added to the radial coordinate as a drift coefficient here a function, the 
other variables such as the velocity magnitudes in different coordinates are 
all simply zero mean < 9ij >= 0, constant drift < % >= b or linear drift 
terms corresponding to accelerated motion < $ij >= c9ij + b; Note that a 
constant drift value mean is useful in uniform motion descriptions, and a 
linear < dij >= cOij + b drift term corresponds to a constant driving force 
term. 



6 Stochastic Derivation 

The derivation of the PDFs of interacting molecules can be made from 
the stochastics...the deterministic dynamics are the force and potential here 
the Lennard- Jones -like potential. The force after variation via the Euler- 
Lagrange equation is for the radial coordinate and here we mention that the 
different species of molecules are accounted for by the center of mass and 
coordinate separation r = xi — i 2 and M = 1711,712 and 

M dv L = _dLJ(rl 

dt dr ( J 

which can be solved for the radial coordinate and fluctuations related to 
the temperature and thermodynamics added by the diffusion-temperature 
Einstein relation, and since v = ^ 

dr = a(r, t)dt + VDdW(t) (22) 

and here the velocity has been integrated and the drift coefficient is the 

integral w.r.t. time of the differential of the Lennard- Jones -like poten- 

t 

tial a(r,t) = f — 9 ^f^} dt' and the fluctuations are a Wiener process' dif- 
o 

ferential , and is delta correlated as a Gaussian white noise for simplicity 
< W{t)W{t') >= 5(t — t')... future research can generalize the fluctuations 
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to the nonextensive statistics for very good Xi-squared experimental data 
fitting and therefore accuracy. 



The stochastic differential equation Eq.(22) is equivalent to the Fokker- 
Planck equation 

this equation is solved by the short time transition probability 



ID(t-t') 

P(^,M')^ ^ D(t _ v) (24) 

The equation is also solved by the short time transition probability of the un- 

primed arguments a(rv,-,t) drift coefficient . The two-point short time transi- 
ts -a' Vjj.t)) 



tion probability can be set to a one-point probability pfrij, t) = — 
and this is compared to the macroscopic maximum entropy information the- 
ory least biased distribution derived PDF.... the diffusion coefficients being 
different in the square root of the potential, the powers of the //, v as due to 
the partial differentiation of the force term, and the Lagrange multiplier of 
the potential term , which now we see can be obtained by inspection from 
this approach. 



13 



7 Computation and numerical simulation. 

The numerics of computation for simulating chemical catalyzed glycol sys- 
tems are simplified by these derived exact solutions. The distributions are 
then utilized with Monte Carlo utilizing PDFs or stochastic trajectories uti- 
lizing the stochastic differential equations and numerical computed results 
obtained. As the distributions are exact solutions for approximated dynamics 
of interactions of the Lennard- Jones type, results can be calculated directly 
and ensemble simulations performed efficiently Also, the NIST challenge, of 
obtaining a theoretical description for dissolved glycol in water and subse- 
quent thermodynamic and phase properties can be either described analyti- 
cally or simulated by numerics. 

The parameters that are to be obtained either from fitting to experi- 
ment data or by alternative means of energetics theory are the powers of 
the potential energy, //, v of the Lennard- Jones potential LJ(r^) = 4- — %- 

and the Lagrange multiplier that sets the units of proportionality for the 
potential energy and coordinates.... this last is a constant presumably as the 
time dependence is included with the variance or second moment's Lagrange 
multiplier, this can also be determined from experiment or energetics of 
thermodynamical relations. 

8 Conclusion 

In this article we have derived a general //, v form of the 6, 12 Lennard- Jones 
-like fluid theory of Glycol water fluid statistical dynamics from an informa- 
tion theory and equivalently maximum entropy approach. This derivation 
is of the theoretical approach of modeling Hydrogen bonds interactions by 
Lennard- Jones and Lennard- Jones -like potentials. We have derived analo- 
gously to the fully quantum mechanical theory of phase transitions of fluids, 
vapors and solids and the matrix theory of interaction-formation of molecules 
from atomic components a matrix of probability distributions and therefore 
of the observables of the macroscopic dynamics of ether or ethylene and wa- 
ter molecules formed from Hydrogen, Oxygen, Carbon etc. that are in a 
particular fluid interaction phase at any given time. 

This is derived from first principles of atoms forming by interactions the 
molecules of the water, atoms forming the dissolved adducts of the glycol... the 
theory includes vapor-fluid transitions or evaporation. We then solve analyt- 
ically a particular distribution function of the interacting Hydrogen bonds' 
Lennard- Jones -like potential and obtain analytic forms that are normalized. 
This analytic form is due to the radial dependence of the Lennard- Jones po- 

14 



-5 -4 -3 -2 -1 







G 



-1 

-2 
-3 
-4 
-5 
-G 

Fig. 1. p(x,t) vs. x , with y=p(x,t) and x= r»j, the diffusion coefficient 
D = 0.001, and time t = 1 and \i — 6, v — 12, a, b — 1 for a Lennard- Jones 
potential estimation of the distribution behavior... note that as x— > in 
arbitrary units the distribution decays rapidly to zero as expected, and as x 
becomes larger the probability increases, reaches symmetrically a maximum 
, and with increasing separation the probability of any two molecules 
interacting again decays to zero, though slowly with separation distance. 
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tential and we compare the macroscopic derivation to the stochastic deriva- 
tion by which the equivalent PDF with drift coefficient is obtained. We then 
discuss future research directions of computational simulations performed on 
fluids which the dissolved glycol water fluid is expected to have reduced freez- 
ing point, and which has been reported elsewhere and as cited here [13] which 
we will report in detail on for materials of water and glycol for simplicity. 
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